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ABSTRACT 

I analyse the maps recording the travel-time shifts caused by averaged plasma anomalies under an 
"average supergranule" , constructed by means of statistical averaging over 5582 individual supergranules 
with large divergence signals detected in two months of HMI Dopplergrams. By utilising a three- 
dimensional validated time-distance inversion code, I measure the peak vertical velocity of 117±2 ms _1 
in depths around 1.2 Mm in the centre of the supergranule and root-mean-square vertical velocity of 
21 ms _1 over the area of the supergranule. A discrepancy between this measurement and the measured 

surface vertical velocity (a few ms _1 ) can be explained by the existence of the large-amplitu de vertical 

flow under the s urface of supergranules with large divergence signals, recently suggested by |Duvall fc | 



Hanasoge (20121. 



Subject headings: Sun: helioseismology — Sun: interior 



1. LARGE-MAGNITUDE SUBSURFACE SUPERGRANULAR 

FLOWS? 

The nature of supergranules - convection-like struc- 
tures observed in the solar p hotosphere - is largely stil l 



in debates (see a review by Rieutord & Rincon 2010) 
Only the surface properties ol supergranules are well es 
tablished. The plasma flows within the supergranules are 
predominantly horizontal with root-mean-square velocity 
^300 ms -1 . The vertical component of the supergranu- 
lar flow is difficult to measure, its amplitude is usually 
within measurement error bars. Statistically it has been 
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known. Attempts were made using local helioseismic 
methods with controversial results. Some studies (e.g. 



Duvall 1998 Z hao fc Kosovichev||2003| to name a few) re 
vealed supergranules as convection- like cells ext ending to 
depth s of 8 -25 Mm with a de e p "ret urn flow" . JWoodard 
( |2007[ ) and |Jackiewicz et al.] ( 2008[ ) did not detect the 
: low reversal, but pointed out that any inversion for the 
flow snapshot deeper than 4-6 Mm is dominated by the 
random noise and thus does not re veal any in f ormat ion 
about the deep supergranular flow. Hathaway (2012) re- 
cently showed that supergranules may extend to depths 
equal to their widths indicating that analysis of super- 
granules with differing sizes may lead to a different depth 
structure. Numerical simulations of near-surface Sun-like 



convection (e.g. Ustyugov| |2008 ) indicated that the am- 
plitude of the flows on supergranular scales decreased 
with increasing depth, although the flows were highly 
structured on smaller scales. 
Recently, a ve ry surprising result came to light from 



helioseismology - Duvall & Hanasogel (|2012[) claimed that 



the careful analysis of travel-time maps provided evidence 
for a large-amplitude flow under the surface of supergran- 
ules. No inversion was performed in this work. The 
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authors detected a systematic offset in the travel-time 
shifts measured for large separations between the mea- 
surement points using a special spatio-temporal filtering 
of the data. This non-standard procedure was selected 
to avoid the cross-talk between the horizontal and verti- 
cal flow in supergranules, which was known to be hard 
to avoid in inverse methods. They concluded that the 
offset was due to the large-amplitude vertical flow below 
the surface. Using a simple Gaussian model, constrained 
by surface measurements, they predicted a peak in the 
vertical flow of 240 ms _1 at a depth of 2.3 Mm. 

In this study, I use a recently implemented, improved 
time-distance inversion code (which minimises the cross- 
talk between the flow components, the main issue leading 
|Duvall iE~ Hanasoge] to use non-standard time-distance 
methods) to verify their hypothesis. 

2. INVERSION SETUP 
The code used in this study (described and tested in 



Svanda et al. 2011) implem ents Subtractive Optimally 
Localised Averaging (SOLA; Pijpers fc Thom pson 1992| 
three-dimensional time-distance ( |lju^nreT al.||l993[ ) m- 
version. The SOLA algorithm allows to balance between 
various terms in the cost function. Among those, im- 
portant terms evaluate the localisation in the Sun (an 
averaging kernel), the level of random noise in the re- 
sults, and cross-talk contributions. By choosing values 
of the trade-off parameters, one can regularise strongly 
about some of these terms, while relaxing demands on 
the others. 

The properties of the hypothetical large-amplitude 
supergranular flow require a special setup of the time- 
distance inversion. To avoid misleading signals from the 
side-lobes of the averaging kernels, one has to pick the 
solution with good localisation and therefore higher noise 
level. Luckily, one can use the statistical approach - the 
trade-off parameter is setup so that the averaging kernel 
is highly localised, while the noise level is too large (of the 
order of 100 ms^ 1 ) to do the tomography of individual 
supergranules. By averaging N supergranules (containing 
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FIG. 1 Two examples of average difference travel- 
time maps used in the inversion. The /-mode 
travel times with centre-to-annulus and centre-to- 
quadrant geometry with annulus radius 14.6 Mm 
are presented. The out-in measurement is sensi- 
tive to the vertical flow and the divergence of the 
horizontal flow, while west— east measurement is 
practically sensitive to the west-east component 
(v x ) of the horizontal flow only. 

independent realisation of the travel-time noise) the noise 

level is reduced VW-times. 

The SOLA inversion was performed in a Cartesian ap- 
proximation (x and y for the horizontal coordinates and 
z for the vertical) using ridge-filtered sensitivity kernels 
(modes / to p%) with centre-to- annulus and centre-to- 
quad rant averaging geometries ( |Kosovichev k, Duvall 
19971. The radii of the annuli ranged from 5 to 2U pixels, 
with a pixel size of 1.46 Mm (0.12°). The ridge filtered 
sensitivit y kernels were compu ted using the Born approx- 
imation ( |Birch fc Gizon|2007| ). Only the near sub-surface 
layers were targeted. 

2.1. Travel-times 



In the current study, I used data similar to Duvall fc] 
Hanasoge (2012). Dop plergrams from Helio seismic and 
Magn etic Imager (HMI; |Scherrer et~aL|2012[ |Schou et al" 
2012) were used to produce 64 consecutive twelve-hour 
datacubes for the travel-time measurements, starting on 
10 June 2011. The inversion was designed to be applied 
to travel ti mes measured from Michelson Doppler Im- 
ager (MDI; Scherrer et al.fl995" ) full-disc Dopplergrams. 
Therefore, the HMI Dopplergrams were degraded to MDI 
full-disc resolution by means of down-sampling and appli- 
cation of a wave-vector filter, so that the power spectrum 
of the degraded HMI datacubes was close to the power 
spectrum of MDI full-disc datacubes. Only the disc- 
centre region was tracked in each datacube; the tracking 
and mapping was done using the code drms.tracking 
(Schunker & Burston, unpublished) implemented in Ger- 
man Science Center for SDO. The Dopplergram datacubes 
were spatio-temporally filtered with ridge filters, separat- 
ing /, pi, and P2 modes. The tr avel times were measured 
following Gizon fc Birch| (2004), using centre-to-annulus 
and centre-to-quadrant geometries with radii of the annuli 
5 to 20 pixels (consistently with the sensitivity kernels 
used in the inverse problem). 

3. AN AVERAGE SUPERGRANULE 

To reduce the level of random noise in the flow maps 
retrieved by good-localisation inversion, I constructed a 
set of travel-time maps representing an 'average super- 
granule', thus reducing the level of random noise in travel 
times. 

3.1. Identification of supergranules 

The locations of individual supergranules w ere deter- 
mined using an al gorithm similar to that used by |DuvalQg| 
Birch| ( |2010[ ) and |Duvall fc Hanasoge| ( |2012[) . Supergran- 
uies were searched tor in the centre-to-annulus difference 
travel-time maps measured from /-mode filtered dat- 
acubes for the distance range of 8-12 pixels (0.96-1.44°). 
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The travel-time maps were smoothed with a Gaussian 
window with Full-Width-at-Half-Maximum of 3 Mm to 
suppress random fluctuations. Such maps were sensitive 
to the horizontal divergence of the flow assumed to occur 
inside supergranular cells. 

1. All pixels, where the sound travel time was smaller 
than the mean by more than two standard devia- 
tions, were tagged. This marked the regions with 
large divergence signal, where the locations of the 
supergranular cell centres could be expected. 

2. For each tagged point a square neighbourhood (with 
a side length of 30 Mm) was taken and the position 
of the maximal divergence (thus minimal travel 
time) in this small region was determined. The 
location of this pixel was tagged. As a consequence, 
several points were possibly tagged in the large- 
divergence region of the same supergranule. 

3. All tagged points were cycled and the distance to 
all other tagged points was calculated. Should the 
distance between any pair of points be less than 
23 Mm, only the point with the minimal value of 
the travel-time was preserved. Only one point then 
represented the original cluster of points. This step 
removed duplicities. 

In all 64 datacubes, the procedure identified 5582 super- 
granule centres. 

3.2. Near-subsurface flows 

For the inverse modelling, the set of all 64 travel-time 
maps was averaged about the location of the supergran- 
ular centres. This resulted in a single set of travel-time 
maps, representing travel-time shifts caused by averaged 
anomalies under the average supergranule. Examples are 
displayed in Fig. [T] 

The averaged travel-time maps were independently in- 
verted for all flow components. The averaging of travel- 
time maps over 5582 supergranules reduced the noise 
level by a factor of \/5582 = 76, thus having the level of 
random noise in the inverted velocities less than 2 ms _1 . 
Averaging kernels for the vertical component of the flow 
are displayed in Fig. [5] (upper panel). The flow estimates 
are averaged mostly over depths of 0-2 Mm with some 
contribution from deeper layers. The centre of gravity of 
the averaging kernels is located at a depth of 1.2 Mm. 

The resulting flow field in the average supergranule at 
depths of around 1.2 Mm is displayed in Fig. [3j While 
the horizontal components of the flow velocity agree with 
what is expected from the literature, the magni tude of 
the vertical flow is very large, in agreement with Duvall 



Inversions for Average Supergranular Flows 
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FIG. 2 Vertical cuts through averaging kernels for 
the good-localisation (upper) and low-noise (lower) 
inversion of near-subsurface vertical flows. Hori- 
zontally, the averaging kernels are roundish. Over- 
plotted contours, which are also marked on the 
colour bar for reference, denote the following: half- 
maximum of the kernel (white) and 5% of the 
maximum value of the kernel (grey). Important 
component of the averaging kernel {K\) is better 
localised in near sub-surface layers in the case of 
good-localisation inversion. Noise levels of these 
inversions are discussed in the main text. It is also 
evident that the contribution of the cross-talk is 
negligible as the cross-talk averaging kernels {K.% 
and fCy) are practically zeroes. 
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& Hanasoge ( 2012 1 . Note that the inverted flow roughly 
lulhlls mass conservation. 

The magnitude of the vertical and horizontal flow, as a 
function of the distance R from the centre of the average 
supergranule, is plotted in Fig. [4| The vertical velocity 
peaks in the cell centre with a value of 117 ms -1 . The 
root-mean-square vertical velocity over the area of the 
average sup ergranule is 21 m s -1 , well within the mea- 



surement by |Hathaway et aL (2002). The second peak in 
the vertical velocity, at a distance of 38 Mm from the cell 
centre, indicates the average location of the neighbouring 
supergranules in my sample. 

The known location of the supergranules allows one to 
also average the measured Dopplergrams and obtain the 
average line-of-sight velocity, which corresponds to the 
vertical yelocitv at the surface. My results are similar to 



those of |Duvall fc Birch| ( |2010[ ), i.e. 12±1 ms -1 in the 
cell centre. 

3.3. Comparison to regular sub-surface flow field 

To check the consistency of my inverse modelling, I 
ran another inversion with the same setup and the same 
target depth, only regularising more strongly about the 
random-noise level term. This approach allows one to 
"see" a snapshot of the near sub-surface flow field on 
supergranular scales. An example flow map achieved 

J100 




Fig. 3. — The flow field in the average supergranule at depths 
around 1.2 Mm. Vertical component (v z ) is indicated by colours, 
horizontal component by arrows. The largest arrow corresponds to 
the flow of 250 m s — 1 . The large annulus of positive vertical velocity 
denotes the average location of the neighbouring supergranules. A 
dashed line indicates the half-way to the neighbouring supergranules, 
thus indicating the size of the average supergranule. 
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with this approach is displayed in Fig. [5j Obviously, no 
lOO-ms -1 vertical flows are present in such a map with 
vertical velocity peaking at around 10 ms -1 , well within 
the limits established by previous studies. Regularisa- 
tion about the noise term in the inversion cost function 
leads to the poorer localisation in the Sun, in particu- 
lar leading to wide lobes in the averaging kernel (see 
Fig. [2]). Application of the low-noise inversion to the 
average-supergranule travel-time maps provides a peak 
vertical velocity of 5 ms^ 1 , twenty-time less than the 
application of the good-localisation inversion. 

The obvious discrepancy can be explained by the dif- 
ferences in averaging kernels. To verify this hypothesis, 
I construct a simpl e model of the supergr anul ar vertical 



flow, following bo th Duvall & Birch 
Hanasoge (2012). The supergranu 
(R, z) is modelled as 

v z (R,z) = vnJo ( 27r ^J ex P 



p010| ) and [Duvall fc 
ar vertical flow field 



2(7, 



(1) 



where the horizontal part is approximated by a Bessel 
function of the first kind with L = 15.1 Mm being the 
half-si ze of the supergranule adopted from |Duvall fc Birch| 
(20101, and the overall magnitude is scaled by a Gaussian 



vertical envelope with vq — 240 ms , zo = —2.3 Mm, 
an d a z = 0.912 Mm, followin g the simple model used 
by Duvall & Hanasoge (2012). In this approximation, 
the vertical flow peaks at the value of 240 ms -1 at the 
depth of 2.3 Mm and decreases towards the surface, with 
a magnitude of 10 ms -1 at surface levels. 

I convolve the model for the vertical flow with the av- 
eraging kernels resulting from both the good-localisation 
and low-noise inversions. The resulting maps do not con- 
tain any realisation of the random noise, but, due to the 
expected noise levels less than 2 ms -1 , it can be con- 
sidered negligible. The resulting peak vertical velocities 
that are retrieved by the inversions from the modelled 
flow field are 26 ms -1 in for case of the good-localisation 
inversion, and 2.8 ms -1 for the case of the low-noise 
inversion. 

The values retrieved from the model do not agree with 
the values retrieved from the travel-times inversion. There 
is one order of magnitude difference which can sufficiently 
be explained by the different averaging kernels. This order 
of magnitude difference agrees well with the order of mag- 
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Fig. 4. — Azimuthal average of the radial (with respect to the 
cell centre, red diamonds) and vertical (v z , blue squares) near- 
subsurface velocity in the average supergranule. The line-of-sight 
velocity [vi os with negative sign; green full circles) obtained from the 
stacked Dopplergrams of supergranules is over-plotted for re ference. 
The green curv e is practically identical to that presented in |Duvall| 
Birch| (|2010|) . Error bars for all plots are 2 ms -1 or less. Results 



ram otner studi es are also overplo tted for reference: the best-fit 
model of IDuvall &: Hanasogel (|2012| at the depth of 1.2 Mm is with 
the black X and the maximum range of the vertical flo ws measured 
in supergra nular centres at the depth of ~ 1 Mm by Jackiewicz 
et al. ( 2008 ) with the black open circle and error bar. Note that at 
the depth of 1 Mm the best-fit model of |Duvall fc Hanasoge| l |2012| 
predicts the vertical flow of 76 ms^" 1 . 
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Fig. 5. — An example of the near-subsurface flow field obtained 
by the low-noise inversion. The measurement were averaged over 
1 day, the predicted errors are 34 ms _1 in horizontal components 
and 5 ms _1 in vertical component. The scale arrow denotes the 
horizontal flow of 250 ms _1 . 

nitude difference obtained from the average-supergranule 
travel-time maps. This naive experiment demonstrates 
that the existing large- amplitude flow might be hidden 
for some inversions with extended averaging kernels. 

4. CONCLUDING REMARKS 

The very large amplitude of the vertical flow measured 
by time-distance helioseismology is somewhat surprising. 
Not indicating possible problems in the forward and in- 
verse modelling, there are three main issues that could 
lead to biassed results. 



• Identification of supergranules The method 
used here prefers the supergranules with highest 
divergence signals. It detected only 5582, while 
something like 50000 is expected, when approximat- 
ing the supergranule by the cell with characteristic 
size 30 Mm in a hexagonal pattern covering the 
whole Sun. Thus, the use of a better segmenta- 
tion algorithm to identify all cells may bring a new 
insight into the problem. 



• Stacking of supergranules The detected super- 
granules are stacked to the point of the maximal 
divergence, which does not necessary need to be in 
the geometrical centre. To overcome this issue, I 
modified the detection code to fit paraboloids in the 
vicinity of the points with the highest divergence 
signals and to choose vertices of these paraboloids 
as locations of the cell centres. This did not change 
the results at all. 



• Foreshortening In this study, all safely detected 
supergranules lying within 24 degrees from the disc 
centre were used to stack the average supergranule. 
In the sensitivity kernels, foreshortening is not ac- 
counted for. Limiting the location of supergranules 
used of 15 degrees in heliocentric angle, only 1014 
supergranules were detected. The amplitude of the 
vertical velocity did not change at all. 



I conclude that there are supergranules on the Sun 
(around 10% of the expected count of all supergranules) 
which exhibit large divergence signals connected with 
the large-amplitude subsurface flows. These flows are 
probably very narrow and might have been missed in 
the past by inversions with a low noise level preference. 
The model of the "average supergranule" inferred here is 
representative only to this small fraction of all expected 
supergranules and does not have to describe the true rep- 
resentative cell. To do so, a better segmentation method 
must be used to construct a model of a "true average 
supergranule" . 
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